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In this second paper, we develop an analytical theory of quasi-equatorial lensing by Kerr black 
holes. In this setting we solve perturbatively our general lens equation with displacement given in 
Paper I, going beyond weak-deflection Kerr lensing to third order in our expansion parameter e, 
which is the ratio of the angular gravitational radius to the angular Einstein radius. We obtain new 
formulas and results for the bending angle, image positions, image magnifications, total unsigned 
magnification, and centroid, all to third order in e and including the displacement. New results on 
the time delay between images are also given to second order in e, again including displacement. 
For all lensing observables we show that the displacement begins to appear only at second order in 
e. When there is no spin, we obtain new results on the lensing observables for Schwarzschild lensing 
with displacement. 
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I. INTRODUCTION 

In Paper I jl| we derived a general lens equation and magnification formula governing gravitational lensing by Kerr 
black holes. Our equation took into account the displacement that arises when the light ray's tangent lines at the 
source and observer do not intersect on the lens plane. 

In this second paper we study the new lens equation of Paper I. We shall restrict our attention to quasi-equatorial 
lensing by a Kerr black hole, and address lensing observables in this regime. Our Paper II begins by obtaining 
the full light bending angle with "horizontal" and "vertical" components for an equatorial observer and light rays 
that are quasi-equatorial. Next, we develop analytically quasi-equatorial Kerr lensing beyond weak-deflection Kerr 
lensing to third order in e, which is the ratio of the angular gravitational radius to the angular Einstein radius. 
Specifically, we solve our lens equation perturbatively to obtain formulas for the lensing observables: image position, 
image magnification, total unsigned magnification, centroid, and time delay. It is shown that the displacement begins 
to affect the lensing observables only at second order in e, and so can safely be ignored for studies of first-order 
corrections to weak-deflection quasi-equatorial Kerr lensing. Finally, the findings in the paper also yield new results 
on the lensing observables in Schwarzschild lensing with displacement. 

The outline of the paper is as follows. In Section|lT]we carefully list the assumptions used throughout the paper. In 
Section [Till we determine the "horizontal" component of the bending angle; the details of the computation itself appear 
in Appendix |B] Section IIVI contains the bulk of our results: we derive formulas for image position, magnification, 
total magnification, centroid, and time delay, all expressed perturbatively to second or third order in the expansion 
parameter e. Due to the lengthy forms of the third-order results, we list them in Appendix |D] Also, as a simple 
application of our results, we plot the image correction terms as a function of the angular source position and compare 
with the Schwarzschild case. Appendix IB 31 contains new results about the vertical bending angle in quasi-equatorial 
Kerr lensing. 
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FIG. 1: Geometry for lensing with displacement by a Kerr black hole. Taken from Fig. 3 of Paper I [ll]. 



II. DEFINITIONS AND ASSUMPTIONS 

We define Cartesian coordinates {x,y,z) centered on the Kerr black hole and oriented such that the observer Ues 
on the positive a;-axis. As in Paper I, we introduce spherical polar angles centered on the observer and defined with 
respect to the optical axis (the x-axis) and the lens and light source planes. The vector to the image position is 
described by the angle -d it makes with the optical axis and by an azimuthal angle if in the lens plane. Recall that 
■i? is strictly positive and within the interval (0,7r/2), while < if < 2tt. Similarly, the vector to the source position 
is described by the angle B it makes with the optical axis and by an azimuthal angle x in the light source plane; B 
is also strictly nonnegative and within the interval [0,7r/2), while < x < 27r (see Fig. 1 of Paper I). A diagram of 
the lensing geometry is shown in Fig. [TJ Note that — 7r/2 < -dg < 7t/2 and < ips < 27r. We adopt the same sign 
convention for 'ds as in Paper I: if 1^5 is measured toward the optical axis, then it will be positive; otherwise it is 
negative (e.g., the -ds shown in Fig. [1] is positive). 

We now state the assumptions under which we will work throughout this paper: 

1. The Kerr black hole and the light source are not at cosmological distances, so that ds = cIl + dLS, where ds 
and are the perpendicular distances from the observer to the source and lens planes, respectively, and c?ls 
is the perpendicular distance from the lens plane to the source plane; 

2. Both the source and observer are in the asymptotically flat region of the Kerr spacetime, and the observer lies 
in the equatorial plane of the Kerr black hole. This last condition implies that the coordinates (x, y, z) coincide 
with the Boyer-Lindquist coordinates {X,Y,Z) centered on the black hole (see Fig. [3] below); 

3. The source is not required to be incrementally close to the optical axis and can be either on the equatorial plane 
or slightly off it, so that x ~ Xo + <5x, where xo = or tt. Similarly, the lift of the light ray off the equatorial 
plane is small, so that = ifQ + 6ip and ips = tpo + + ^^S: where ipo — (retrograde motion) or tt (prograde 
motion), and where Stp and Stps are small and considered only to linear order. We henceforth refer to this as 
the quasi-equatorial regime. 



In Paper I we derived the following general lens equation governing lensing by a thin deflector, for source and 
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observer in the asymptotically flat region: 



ds tan B cos X = dLtawd cos^p + d^stani^s cos(p5 + dj, , (1) 
dstanSsinx = ditanz^sin^j + dLstands sin ips + . (2) 

Here the displacements are shown explicitly; note that (d^ + d^) = d in Fig.[T] 

Specializing to the case of an equatorial observer in the Kerr metric, we also derived in Paper I the following lens 
equation with the displacements implicitly included: 



^5 tan S cos X 
ds tan B sin x 



d stands cos if s 
dLstan-ds sinips 



sin "d cos ip 
cos -ds 
dL sin ■& 



1 — sin i?5 sin ips 



(3) 
(4) 



cos (f sin ds tan ds sin ips cos ips + (sin^ ip — sin^ ds sin^ (fs) 



1/2 



This is valid for all light rays, whether they loop around the black hole or not, as long as they lie outside the region of 
photon capture. No small-angle approximation is required. Comparing our two lens equations allowed us to extract 
the displacement terms: 



dy = dL sind cos ip 
dz = —dLtandsinip 



1 



1 



cos ds cos d ^ 
c?L sini? 



1 — sin^ ds sin^ ips 

1 /2 

COS (p sin 13 s tan i^s sin (ps cos (ps + (sin^ (^9 — sin^ i^s sin^ (fs) 



(5) 



(6) 



III. QUASI-EQUATORIAL KERR LIGHT BENDING 



With that as background, we begin this paper by calculating the component of the bending angle in the equatorial 
plane, which is the xy-plane in Fig. [T] we call this the "horizontal" component. Due to the technical nature of the 
calculations, we quote the key results here and refer to Appendix IB] for the detailed treatment. Note from Fig. [1] that 
according to the way the angles d and ds are defined, they may lift off the xy-plane or a plane parallel to it. Let us 
define d and ds to be their projections onto the xy-plane, respectively. Without loss of generality, we choose the same 
sign conventions for d and ds as we chose for d and ds, in which case we can unambiguously write the "horizontal" 
component of the bending angle as 

Cihor ^d + ds . (7) 
Note that the positivity of d and the fact that the bending angle is nonnegative forces the condition 

ds>-d. 

(Indeed, with our signs conventions the condition ds < ~d would be equivalent to repulsion of the light ray.) Writing 
d and ds in terms of the angles d, ip, ds,ips, we have 

d — tan""'^ (tan?? cos 1^9) , (8) 
ds — tan~"^(tani?5 cos(7r — (^5)) . (9) 

As stated in assumption (3) above, in the quasi-equatorial regime we have 

(p = ipQ + S(p, ips ^ (po +Tr + Sips , (10) 

where ipo is either (retrograde motion) or tt (prograde motion), while Sip and Sps are small and considered only to 
linear order. In this regime eqns. (|S]) and ^ simplify to 

d K ±d , 
ds ~ ±ds ■ 
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Since d and i?5 have the same signs as i? and ^s, respectively, we discard the negative solutions, so that eqn. ([7]) 
reduces to 

dhor«^? + l?S, l9s>-^- (11) 

Thus in the quasi-equatorial regime we may use the full angles and "ds in place of their respective projections onto 
the ccy-plane. With that said, we show in eqn. (jB19|) of Appendix [B] that the "horizontal" bending of light has the 
following invariant series expansion: 

a,Ub) - A,(^) + + A,(^Y + a4^Y + 0(^)\ (12) 

where 



b / ^\ b J "V6/ \ b J \ b 



Ai = A, (13) 
IStt 

A2 = — -4sa, (14) 
128 

As = — - 107rsa + 4a2 ^ (15) 

O 

34657r _ ^ 2857ra2 ^ „ 

^4 = — 192sa+ 4sa^ . 16) 

64 16 

The variable s equals ±1 depending on whether the light ray undergoes prograde (+1) or retrograde (—1) motion (see 
eqn. (jB5|) in Appendix [B] below) . Note that eqns. (fT2 |) - (fT6l) are consistent with the bending angle obtained by Iyer 
& Hansen [l5| by a different means — note also that their bending angle is consistent with the exact bending angle 
[l6j . We remind the reader of our conventions in Paper I. The parameter m, is the gravitational radius and a is 
the angular momentum per unit mass, 

GM, J, ,^„, 

= — 5— , a = — — , (17) 
c2 cM, 

where A^, is the mass of the black hole and J^, its spin angular momentum (see, e.g., [l^, pp. 322-324]). Note that 
both m, and a have dimensions of length. The quantity a is a dimensionless spin parameter: 

a 

m. 

Unless stated to the contrary, the black hole's spin is subcritical; i.e., a? < 1. Finally, b = d^sini? is the impact 
parameter (see eqn. (A8) in Paper I), which is a constant of the motion. 

When there is no spin, the coefficients reduce to Ai = A, A2 = 157r/4, A^ = 128/3, and A4 — 34657r/64 and recover 
the Schwarzschild bending angle in [13. Also, eqn. ([T^ shows that in the weak-deflection limit (at first order in m,/6) 
the Kerr bending angle agrees with the Schwarzschild bending angle. The spin enters only in higher-order correction 
terms. The sign is such that the spin makes the bending angle larger for light rays that follow retrograde motion 
(s = —1). This makes sense intuitively because retrograde rays spend more time in the presence of the black hole's 
gravitational pull. 



IV. OBSERVABLE PROPERTIES OF LENSED IMAGES 



In this section we derive asymptotic formulas for image position, image magnification, total unsigned magnification, 
centroid, and time delay for quasi-equatorial Kerr lensing with displacement. 



A. Quasi-Equatorial Lens Equation 

We begin with our general lens equation p^ - ([T^ and insert a bookkeeping parameter ^ to monitor the displacement 
in either dj, or d^: 

ds ta.nB cos X = dL tan 1) cos (p + c^ls tan i?5 cos 1/55 -|- ^ dy , (18) 
dstaniJsinx = tani? siniy9 + tani^s sint/jg -I- ^ . (19) 
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(The displacements dy and are given by eqns. ^ and ([6]).) We can take ^ = 1 to include the displacements 
properly, or choose ^ = if we wish to ignore the displacements (in order to connect with work in [13, [isj). 

Beginning with eqn. (fT8| . we substitute eqn. ([5|) in place of dy and Taylor expand in the small angles Sip, S(ps: and 
Sx, to obtain 



qtsniB = {1 - D)tan'd - Dta.n'ds + ^(1-D)sin?9 



1 



1 



cos COS i9 



0{2) 



(20) 



where q — cos(xo ~ fo), D = d^s/ds-, and O (2) indicates terms that are second order in Sip, Sips, and/or Sx- (Below, 
we incorporate the sign q into the tangent so that the left-hand side is written as tan(qi3) and we think of qB as the 
signed source position.) This is the "horizontal" component of the lens equation. Bear in mind that ^ identifies terms 
associated with the displacement. Including the displacement by setting ^ = 1 in eqn. (^0)) yields 



q tanS = (1 - £>) 



sin?? 



Dtan'ds + O (2) 



cosds 

Thus, to lowest order in out-of-plane motion we recover the same lens equation as in the Schwarzschild case (see 
eqn. (19) of Paper I). 

We use ds = cihor — taking dhor from eqn. ((T2|) . and introduce scaled angular variables: 



qB 



m. 

dL 



t?, = tan"^ 

Here the natural angular scale is given by the angular Einstein ring radius 

^E = 



£ = 



AD 



AGM.dLs 



Am,D 



(21) 



(22) 



Note that we have defined the scaled source position /3 to be a signed quantity, with a sign that indicates whether 
the source is on the same or opposite side of the lens as the image. In eqn. ([T2|) we wrote the bending angle dhor as a 
series expansion in m,/6. For analyzing the observable image positions, e is the more natural expansion parameter. 
To convert Shor into a series expansion in e, note that according to eqns. ([2T|) and (|22|) and the fact that b = dLsini}, 
we have 



m. 



1 



22AD^ 6*3 
45" 



-0{ef 



(23) 



6 sin(4L>e6i) 9 3 

As in [l7l - [T9| . we postulate that the solution of the "horizontal" lens equation (|20|) can be written as a series 
expansion of the form 



= Co 



^2 



^3£" 



0{ef. 



(24) 



Converting now to our scaled angular variables ((2T|) -p4 | . our quasi-equatorial "horizontal" lens equation ((20)) takes 
the form 





1 




^0 L 



s a - — + (1 + 02), 



(25) 



24 9l 



768 9^ 



12 s a (Stt -A9i)- 24 d^ - 384 + 3 Oi (IStt -Wi) + 24 9^ 6*2 (1 + 6*5) 



6*0 (48 D9a + 9l{-2^^ - 79a + 29^^)) + 19291 (1 - D)(l - 2 Z? 9^) ^ 



128 s d {384 + 96ii(-57r + 29i) + 4:9q (-392 + 89o{2D{-3 + 209^) 



+ 3(1 - D){~1 + £161^) 0) } - 768 d^ {2tt - 301) + 768 s d^ 
+ 157r {-3(487 + mf) + 32 9a (302 + 86'o(2£'(3 - 2D9l) 
+ 3(-l + D){-1 + D9l) 0) } + 256 {30? + 30^(1 + 0^)03 

+ 201 (72 + 0o(-302 + 40o(7:'(-6 - 7i:>0^ + 6i:>0^) + 3(-l + i:>)(l + 2D9l) 0)) } 



+ 0{ef 
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Note that displacement terms (indicated by ^) only begin to appear at second order. Also, since we are simultaneously 
expanding tanqS — tan(4/3/?e), note the occurrence of (5^ in the £^ term. 

Now we turn to the "vertical" component of the lens equation, namely, eqn. (fT9|l . Substituting eqn. (jG]) in place of 

and Taylor expanding in the small angles Sip, Stps, and Sx, we obtain 

((5x)(qtan6) = (Sip) {1 - - D)tSiiiiS - (Sips) D tands (26) 
+ ,^(1 - D) sin?9 |(^</J5 siniJs tani?^ + [{5ipf - (Sips)^ sin^ ^s] ^^^| 
+ 0{2). 

Next, we use eqn. (jB24[) in Appendix [B] to write {S(ps) s'm'ds = (Sf) W{d): 

(5x)(qtanS) = 6ip{l - - D)tan^ + dip DW{^){cosi}s)~'^ (27) 
+ C{1 - D)sin^{dipW{^)taii^s + Sip[l -WiT)f]^/^} . 

(In eqn. (|B26p of Appendix [B] we show that 1 — VF(i?)^ > 0, so this equation is never complex- valued.) Finally, we 
convert to our scaled angular variables ([2T1) -(|24 |) and expand in e, obtaining 



(3 + ^D^/3'e' + 0{ef]Sx^5ip 



(92 



2sa-^-^ + {l + 0l)i 



(28) 



1 



24 



s a (907r - 966li) + 2a^ + 45716*1 - 24(16 + Of) + 86*0 {8D^el{-7 + 26^ + 6^ 



+ 3(02 + 0102 + S0oO - 24D6o{-2 + C + 20^^)) 

aVi + 72a^(-957r + 96 6ii)-|-3072sa^ + V2 



768 6(4 



where 



Vi = 64s [960 + 96'i(-157r + 86ii) + 166io (-3612 + 46io(i:'(-6 + 5i:'6'2)-3(-l + D)(-l + i:i6|2)e))] , 
V2 = IStt [-3(487 + 4861^) + 326io (36I2 + 86lo(2D(3 - 2Del) + 3(-l + D){-1 + D6l)S,))\ 

+ 256 [361? + 30^(1 + 6'2)6i3 + 26ii (72 + 6io(-36i2 + 46io(i:'(-6 - IDOl + QDOl) 

+ i{-l + D){l + 2Del)0))] ■ 

This is the "vertical" component of the lens equation. We will use it to obtain a relation between the small angles 5x 
and 5ip. To that end, we divide eqn. (1^ by eqn. (1^ to eliminate (3: 



5x = 5^ {l + 



a 



a 



7686*3 (6i2 „ 

768sa2(l-36'g + 36l^) + sK 



s (-57r + 46l^(57r - 126li) + I66I1) + 16a(l - 261^) 

24 a (7r(101 - 262 6^ + 221 ^j^) - 64(1 - 3 6i^ + 4 9^) 61] 



(29) 



where 

K = 2257r2(-l + 46*^) - 4807r(l- 36*2 + 86'*)6'i+ 7686*2 

+ 256 [0n {-302 + 0oi-8D^ ^o(-l + ^o)(-5 + 6C) + 24i?(-l + el){-2 + ^ + 29l^) 
- 3(16 + Wl - 8^ + 0o{-'^02 + eo(-16 - ml + 3^2 + 80))))] • 

Observe that in general 5x 7^ Sip in the regime of quasi-equatorial lensing. Thus when a 7^ 0, the light ray's trajectory 
cannot lie in a plane other than the equatorial plane (in which case 5ip — 5x — (^)- Notice that displacement terms 
(indicated by appear only at third order in e. 
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B. Image Positions 



We now solve our "horizontal" lens equation (1251) term by term to find 6*0,^1, 6*2 , and 63. The zeroth-order term is 
the familiar weak-deflection lens equation for the Schwarzschild metric, 



13 ^Oo 



1 



which yields the weak-deflection image position 



(30) 



(31) 



We neglect the negative solution because we have explicitly specified that angles describing image positions are 
positive. For a source with /3 > 0, the negative-parity image is obtained by plugging — /3 in eqn. pip: note that 
eqn. (f3T1) will still be positive. (Note also that we are solving for quasi-equatorial images; there may be additional 
images in the general case.) 

Requiring that the first-order term in eqn. (j25p vanishes yields 



IStt — 16 s a 



(32) 



Likewise with the vanishing of the second-order term, 

1 



64 (6 - D{2D + 6(1 - D)el - DO^)) + 240^ - 12 s a(57r - ^Oi) 



24 00 (^^ + 1) ^ 

-36'i(157r - 86*1) - 192612(1 - D){1 - 2Del) ^ 



(33) 




where we have used eqn. (150]) to substitute for /3 in terms of 6*0. Note that the displacement only affects 62, not 6*0 
and 9i . The third-order image correction term 63 is given in Appendix [D] 

In terms of the source position /?, we can write the terms for the positive- and negative-parity images as 



where we have written to remind ourselves that the two images have different respective values of the pro- 
grade/retrograde sign parameter. In fact, we have = — s+. The terms 9^ and as functions of /3 are similarly 
obtained, but are too lengthy to be written here. Now thinking of the universal relations studied in [l8|, we observe 
that the zeroth-order terms obey 

which are identical to the zeroth-order position relations obeyed by PPN models (see [13 )■ The first-order terms have 

In [l^ it was shown that 9^ + 9^ is independent of source position for static, spherical black holes in all theories of 
gravity that can be expressed in the PPN framework. However, as first shown in (26| . we see that in the presence 
of spin 9^ -\- 9i is no longer independent of source position. This is a direct consequence of the fact that one image 
corresponds to a light ray that follows prograde motion, while the other has retrograde motion. The difference between 
the second-order components is (cf. |18| ) 



9t 



2aV4 + /32 + 
— 307rs+a + a(' 



(16a + 157rs+)v/4 + /32 

8(4 + ^2)3/2 

157rs+) (4 + I 



8(4 + /32)3/2 



256 
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/? (in units of t^s) /5 (in units of &e) 
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FIG. 2: First-, second-, and third- order angular image correction terms as functions of the angular source position /3, for 
a positive-parity image undergoing either prograde (s = -fl) or retrograde (s — —1) motion near the equatorial plane of 
a Kerr black hole. The solid curves represent a Kerr black hole with spin parameter a = 0.7. When a = 0, we recover 
Schwarzschild lensing (dashed curves). For the second- and third- order image corrections, the displacement parameter ^ = 1 
and D = dLs/ds = 0.5. (Note that 61,62, and 63 are dimensionless, but have factors of e,e^, and e^, respectively. Note 
also that (Ils and ds are the perpendicular distances between the lens and source planes and the observer and source plane, 
respectively.) These results hold for a black hole with sufficiently small i9_b. 



where 

I = -4096 + 2257r2 + 20i8D^ + 1607rs+a - 5120^ + 4096D(1 - D)^ 

We can likewise consider 6*^ ± 9^ , but will forgo writing them here. Plots of the image correction terms 6i , 62 and 
63 as functions of the source position /3 are shown in Fig. [51 for a positive-parity image undergoing either prograde 
(s = or retrograde (s = —1) motion. 



C. Magnifications 

In Paper I we derived the following general magnification formula: 



sinS / dB dx dB dx 



sin ■& ydi} dip dip dd 



To compute dB/d-d, we employ the same techniques that led to eqn. (|25p . For dx/d(p, we use eqn. 
dB/dip = for quasi-equatorial lensing.) The result is the following series expansion: 



(34) 
(Note that 



H = fio + ^ii e + H2 £^ + A*3 + O (e) , 



(35) 
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where 



Ml 



M2 



[157r(6l2- 1)2 + 64 s a 6l2]6l3 



16 I 



384(02 _ 1)3(02 + 1) 



+ 1207rs+a(l + IGe^ - 349^ 



1)3 ' 
768 026*^(5 



261 + S^o) 



(36) 
(37) 
(38) 



^ai^o - l)"" (-12288i:i(6io + Oo f + 1024i:)^(l + e'fj-'ief^ + 160'o + 1) 



446'!5 

q3\2 



3961° 



1201") 



2\2^ 



- 36*^ (4096 + 6'2(-6757r2 + 4096(6|2 + 2)))) - 6144(1) - 1) {Oo - 6'g)2(-l + 2i:>0g(2 + 6ig)) ^ 



Note that displacement terms (indicated by ^) begin to appear only at second order in e. The third-order magnification 
term fi^ is given in Appendix |DJ In terms of the source position /?, we can write the terms for the positive- and 
negative-parity images as 



t4 



1 ^ + 2 

2 2|/3|y^ 



IStt/ 



■64s±a 



± 



16 /32(/32+ 4)3/2 ' 

20257r2 - 1024(/32 + 4)(12 4- £1(12 - {(3^ + 18)1))) 



57r s a 
32/32 

16^(1- 



-1± 



384 |/?| (^2 + 4)5/2 
|/3|(/34 + 34;32+48' 



(/32 

D{f3^ 



■f 4)5/2 

+ 6/32 + 6) - 1 



± 2a2 



5^2 + 8 



|/?|3(/32+ 4)5/2 



|/3| (/32 + 4)3/2 



Observe that 



mJ 



A*2 



Mo 

Mo" 

Mr 
Mr 

M2" 



/32 



3|(/32+ 4)1/2 ' 
IStt 

"8(/?2 +4)3/2 ' 



/32(/32 + 4)3/2 ' 

57rs+a 1^1 (^^ -I- 34/32 + 48) 

^6 /?2(^2 _^4)5/2 

57rs+a . 5/32 + 



4a^ 



32ei)(l-i5) , 



32^1-0) 



i:'(/34 + 6/32 + 6)- 1 



16/?2 ' — |^|3(^2 ^„ 4)5/2 

20257r2 - 1024(/32 + 4)(12 -I- £1(12 - [fi'^ + 18)L»)) 
192 |/3|(/32 + 4)5/2 



(^2 + 4)3/2 



(39) 



The expressions /ij^±/i;^ are given in Appendix [P] The zeroth-order sum relation is the same as the universal relation 
found for static, spherical PPN models in [l^. Notice that the zeroth-order difference relation is independent of spin. 
In the first-order difference relation, the right-hand side is zero for PPN models, but nonzero in the presence of spin 
(see also [l^ ) . In the second-order sum relation, the right-hand side is not zero even in the absence of spin. This is a 
consequence of the displacement (indicated by ^). 



10 



D. Critical and Caustic Points 



To determine the set of critical points, we set /i ^ = 0, the reciprocal of the series expansion given by eqn. (j35p in 
Section liV CI and solve for 60,61,62, and 63. This yields the following 0-components: 



9 

^critical 

9 ^ 
1 critical 



1 , (40) 
. 157r 



32 



(41) 



= S-'-^~'-S^ + D-(^4-H]+mH^2)-^, (42) 



-^^2 critical 



32 2048 V 3 



+ /2257r2 „ \ 157ra2 



93f.iticai - -s a — - - 81?^ (1 - - 81? ^ 



V 256 V y g4 

157r(400 - 2257r2 - 4096152(1 - ^) + 2048 ^ - 2048L»(-2 + 3 0) 

8192 



(43) 



where " ± " corresponds to the two values s = ±1. Note that since we are in the regime of quasi-equatorial Kerr 
lensing {tp = ipo + Sip with ipQ — (retrograde motion) or tt (prograde motion)), eqns. (|40p - (|43p do not define a circle 
on the lens plane, but are to be interpreted (by eqn. (IB6p ) as four points (6*^^;^, tt ± Sip), {6~^^^,±Sip) on the lens plane, 
for a given 6(p. We now insert these into the "horizontal" lens equation (l^5t to third order in e and solve for 13. This 
yields the /3-componcnts of the caustic points, which we express here as a series expansion in e to third order: 

a± ^ 57rsa 2 

- ~sas (44) 

16 ^ ' 



caustic 



a 

512 

0{e 



11367ra + s (2257r2 - 4096^ + 4096i:»2(i _ 2^) + 4096i:'(-2 + 3^) 



,4 



The signs ± correspond to prograde (s — +1) and retrograde (s = — 1) motion, respectively. When a — the two 
caustic points converge to one point at the origin of the source plane. Note from the third-order term that the caustic 
points are not symmetric about the vertical axis on the light source plane. 



E. Total Magnification and Centroid 

If the two images are too close together to be resolved (as in microlensing) , the main observables are the total 
unsigned magnification and the magnification-weighted centroid position. Using our results above, we compute the 
total unsigned magnification: 

Mtot = Im+I + = 1^1(^2+4)1/2 - ^2^2'^ 4)3/2 ^ + ~ 1^2) + {^t - ) + 0{Ef, (45) 

where the second- and third- order terms are given by eqns. (|39|) and (|D1I) . respectively. The magnification- weighted 
centroid position (actually, its "horizontal" component, since we are in the regime of quasi-equatorial lensing) is 

_ _ |/3|(/32 + 3) (2-/32)s+a 

iM+l + lMi " /3^+2 + (/32 + 2)2 ^ 

+ 384|m8+X' + /3^)3 -^ + B...,3e3+0(.)^ (46) 

where 

C2,i = 1207rs+a|/3|(2 + /32)(3 + /32)(4 + /32)3/2 + 3g4^2 ^(2-f /32)(-l6-8^2_^/34) 

4(8 + 2/3^ + /?'')] - ^2(2 + [3(6757r2 - 4096(4 + f)) 

-I- 1024(4 + /32) (D(6^2 - £'(-2 + 9/3^ + ^^)) -t- 3(-l + B){-^'^ + 2i:>(6 + 4^^ + ^4)) ^ 

and where the third-order term 0cont,3 is given in Appendix |DJ In 18] it was shown that the first-order corrections 
to the total unsigned magnification and centroid position vanish universally for static, spherical black holes that can 
be described in the PPN framework. In the presence of spin, the first-order corrections are nonzero. Once again, the 
displacement terms (indicated by ^) appear only at second order in e in both eqns. (j45|) and (j46p . 
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F. Time Delay 



In Appendix [Cj we show that the lensing time delay can be written as 

ds 



cr = r(i?,,c) + r(i?obs 



cosB 



where 



-Robs — 



'2 

LS 



. tan^ B) 



1/2 



B = 4/3De , 



and -Robs and Rsrc are the radial coordinates of the observer and source in the Kerr metric. We derive a Taylor series 
expansion for the function T{R) in Appendix [C] (see eqn. (jCll) '). To determine the observable time delay, we evaluate 
T(i?) at i?src and i?obs, and then replace ro with b using eqn. (|B12p . We change to angular variables using h = (11 sini?, 
and then reintroduce the scaled angular variables in eqns. (Pl]) - ([M1) . Finally, we take a formal Taylor series to second 
order in our expansion parameter e. This yields 



T 



l + P'^-el- In 



dL eldl 



Ad 



LS 



IStt — 16 s a 

16^;^ 



1536 6*2 1 



+ 1) 



(47) 



where 



-967rsa(-7 + 



9t) + 768 {291 + 3^0 + ^o) - (1 + ^o) {6757r2 
+ 3072 el [el + 1) (2 + + 6io4 - 2/3^ (6»2 + i) - 4^) 
+ 3072i:»6i2(l + 6i2)(-8 + + 2^25,2 _ 3^4 _^ 4^ ^ g^2^) 
+ 1024£)2(i + ^ (24 - 5/3^)02 + 50^ _ 24^4^) } , 

and the natural lensing time scale is 

dLds 



(48) 



TE 



cdLS 



^E 



Notice that retrograde motion (s = —1) leads to a longer time delay than prograde motion (s = +1), which makes 
sense intuitively. As with our other lensing observables, displacement terms in the time delay (indicated by ^) begin 
to appear only at second order in e. 

The differential time delay between the two images, At = — r+ is such that 



At 

TE 



In 



1536(/32 + 4) 



157r 



+ s+ + 4 



(49) 



where 



D = 967rs+a(4 + /32)(7^2-l)- 768a2|/3|y4T^-|/3|\/4T^{-6757r2(3 + ^2-) 

+ 3072(-8 + 2/32 + 13^) + 10241)2(4 + l3^){18 + 5/3^ - 24^ - 3072D(4 + p'^){6 + - 8^} 



(50) 



V. REMARKS ON LENSING OBSERVABLES 



We make a few remarks regarding our results: 

1. The procedure for solving the lens equations in the quasi-equatorial regime is as follows: given a source whose 
(scaled) location on the source plane is (/3,Xo + ^x)^ we first solve the "horizontal" lens equation (PS)) term by 
term to find 9o,di,92, and 6*3 (all expressed in terms of /3), and then insert these into the "vertical" lens equation 
(P^l and solve for 5(p. The (scaled) locations of the two images in the lens plane are then 

(6*0 + 9fe + 9fe^ + 9fe\ <fo + <5</?=^) , 

where "±" corresponds to s = ±1 and where (^g — for retrograde motion (s = —1) and tt for prograde motion 
(s = +l). 
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2. Note that for all lensing observables — image position, image magnification, total unsigned magnification, cen- 
troid, and time delay — the displacement parameter ^ begins to appear only at second order in e. Therefore 
displacement can safely be ignored for studies of first- order corrections to weak-deflection quasi- equatorial Kerr 
lensing. Note that the displacement affects the caustic positions only at third order in e. 

3. When there is no spin, we obtain new results on the lensing observables due to Schwarzschild lensing with dis- 
placement. Indeed, all of our results in Section |IV] (including the third-order results in Appendix ID]) immediately 
apply to this regime once we set a = and the displacement parameter ^ = 1. This is equivalent to beginning 
with the spherically symmetric lens equation with displacement (given in Bozza & Sereno [3] and eqn. (19) in 
Paper I) and then computing lensing observables perturbatively in e. 

4. If one sets a = = ^ (i.e., if one turns off spin and ignores displacement), then all of our results in Section ITVl 
are consistent with the previous studies of Keeton & Fetters [l3, [3| • 

5. The total magnification and centroid (eqns. (1451) and (j46p . respectively) are consistent with the corresponding 
results obtained in Werner & Fetters |26[ to first order in e. (The analysis in [2^ was carried to first order in 
e and did not consider displacement.) In fact we point out that the "horizontal" and "vertical" components of 
our lens equation (eqns. p5|) and (|28| ) and our magnification terms (eqns. ((36)) and (|37)) ) are all consistent to 
first order in e with those in 26], after an appropriate change of variables. Furthermore, the "horizontal" and 
"vertical" components of our bending angle (see eqns. (jB20l) and (|B28[) in Appendix [B] below) are also consistent 
to second order in e with those in [26| (their bending angles were written to second order in e). 

6. Finally, our image correction and magnification terms are also consistent with those in Sereno & De Luca [2l| to 
first order in e, while the "horizontal" and "vertical" components of our bending angle are consistent to second 
order in s. 



VI. CONCLUSIONS 

In this paper we derived analytical expressions for the lensing observables for the case of quasi-equatorial lensing 
with displacement by a Kerr black hole. We calculated the light bending angle explicitly for an equatorial observer and 
light rays that are quasi-equatorial, and then applied this to a perturbative framework to third order in the invariant 
parameter e, which is the ratio of the angular gravitational radius to the angular Einstein radius. We obtained new 
formulas and results for the fundamental lensing observables: image position, image magnification, total unsigned 
magnification, centroid, all to third order in e, and time delay to second order. Our results made explicit the effect 
of the displacement that occurs when the tangent lines to the rays from the source and observer do not intersect on 
the lens plane. We also showed that displacement effects begin to appear at second order in e, and so can safely be 
ignored for studies of first-order corrections to weak-defiection quasi-equatorial Kerr lensing. Our findings also yield 
new analytical results for Schwarzschild lensing with displacement. 

This analysis has allowed us to go beyond previous work and probe deeper into the gravitational field of a Kerr 
black hole, by providing explicit perturbative analytical formulas showing how each lensing osbervable is affected by 
higher-order terms. Our results should be useful in observing general relativistic corrections, and can also be used as 
a tool in testing Einstein's theory and perhaps also Cosmic Censorship. 
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Appendix A: Transformation from Sky Coordinates to Boyer-Lindquist Coordinates 

In this section we determine the relation between angular coordinates (i?, ip) on the sky as measured by the observer, 
and the slightly modified Boyer-Lindquist coordinates (t, r, C, 0) shown in Fig. [3] 

Recall from Faper I that the latter coincide with the usual Boyer-Lindquist coordinates {t, r, p, (/>), except that the 
polar angle p is shifted to C = 7r/2 — p. To analyze light bending, it is actually convenient to work with another set 
of coordinates, namely, the lens-centered coordinates (r, T, $) shown in Fig. Our goal is to connect the modified 
Boyer-Lindquist coordinate angles {C,4>) to observer-centered angles {djip). This will be done in two stages: first. 



13 




FIG. 3: Cartesian {X,Y,Z) and spherical polar (r, <^) coordinates centered on the black hole, where = 7r/2 — p with p 
the polar angle. The black hole spins about the Z-axis, which corresponds to — n/2, in the direction of increasing (j). The 
equatorial plane of the black hole corresponds to ^ = or the {X, y )-plane. Taken from Fig. 4 in Paper I [ij] . 

by relating (C, </>) to (T,$), and then by relating (T,$) to This simplifies the ray-tracing and the resulting 

geodesic equations. 

Comparing Fig. |3] with Fig. |4] below yields the following relation between {(, (f) and (T, $): 

sin C = sin T sin $ , tan tan T cos $ . (Al) 

In order to relate the observer-centered angles (i?, (p) to the lens-centered angles (T,$), we make the following 
construction. Consider extending the actual light ray to infinity both "behind" the source and "beyond" the observer. 
Note that evaluating at such points is well-defined because the light ray is a linear path in the asymptotically flat 
regions where both the source and observer reside. With that said, the asymptotic "final" angular position the light 
ray reaches is (cf. Section IIIII) 

T/=i?, <^f = ip + TT (prograde motion), 

Ty = — f — Lp (retrograde motion) . 

The asymptotic "initial" angular position from which the light ray originates is 

Tj = TT — "iJs , <^f — Lps (prograde motion), 

T i — IT + 'ds ^ f — 'i^S (retrograde motion). 

Using eqn. (|A1[) . we can find the initial and final positions in terms of the angles {<t),C)- For prograde motion, these 
are: 

sin = sin sin (^5 , sin^j = — sini9sin(^ , (-A-2) 

tan (pi = — tan 'ds cos 1^95 , tan (j)f — — tan d cos . 
For retrograde motion, they are: 

sin = — sin sin (/9s , sin^/ = sini^sint/?, (A3) 

tan (f)i = tan iJg cos ips j tan 4>f = tan cos ip . 
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FIG. 4: Lens-centered coordinates (r, T,<1?). 

We will use eqns. (jA2p and (|A3|) in our derivation of the "vertical" component of the bending angle vector in 
Appendix IB 31 below. 

Appendix B: Quasi-Equatorial Kerr Bending Angle 
1. Equations of Motion for Quasi-Equatorial Null Geodesies 

Recall from Appendix A 1 of Paper I that the equations of motion for null geodesies are 

t = 1 I 2m,r{a^ -sC + r^) 

[a2 +r(r- 2m.)] (r^+a^sin^C) ' 
. _ , [r^ -{Q + t^- z^V + 2m.((£ - zf + Q)r - a^Q] 

"I" — 9,9-2.. ' V^^j 



a2 sin" C 

2am, r + Lr[r — 2m,) sec^ C + tan^ C 



r 



(r- 2m,)](r2 +a2sin2 () 



(B3) 



I ^ C-.WC)V2 ^^^^ 

+ a^ sm (, 

where Q = Q/f ^ and t = C/£, with £ the energy, £ the orbital angular momentum, and Q the Carter constant (a 
and m, are given by eqn. ()17p ) . Now consider an equatorial observer and source in the asymptotically flat region. 
To compute the light bending angle, we focus on null geodesies that remain close to the equatorial plane (which is 
a plane of reflection symmetry) . There are light rays in the equatorial plane that have C = everywhere. There are 
other light rays that remain close to the plane and have |C| ^ 1 everywhere. Such quasi-equatorial light rays must 
have if = ifQ + Sip with ipQ equal to either or tt, and \6ip\ <^ 1. Given the spin configuration, light rays with (po = 
follow retrograde motion and have C < 0, while light rays with tpo = tt follow prograde motion and have C > 0. Thus, 
if we define a sign s by 

_ . _ J +1 prograde motion, (Ti'^\ 

° ^ ^ )—l retrograde motion, ^ ^ 
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then we can identify 

s = — cos (po . (B6) 
We showed in Paper I that the constants of the motion C and Q can be written as 

C = — sin I? cos If , Q. = d\ sin^ d siv? Lp . 
In the quasi-equatorial regime, these become 

C = sb cos 6(p , Q = b^ sin^ Sip . 
We expect 6ip to be of the same order as C, so we can Taylor expand eqns. (IBip - (jB4l) in both ( and Sip. This yields 



t 



( 2 m, a&s 

1 + - 



dr + r(r — 2m, 



0(2), 



± 1 



62 2 m, 52 ^ 

&s / 2m, 

1 



32 ^ _ 2m,) 



0(2), 
0(2), 



C = ±^^{{Sp>f~H'\ 



21 1/2 



0(2), 



(B7) 

(B8) 
(B9) 
(BIO) 



where O (2) indicate terms that are second order in and/or Sip., and we have defined 



F = 1 - s ■ 



1 — s a 



m. 



1 



62 



62 



Notice that r, and do not depend on <^ or Sp at zeroth order or first order. In other words, the "in plane" motion 
is insensitive to small displacements above or below the equatorial plane. By contrast, C lacks a zeroth-order term 
but has a nonzero first-order term. Thus, there is a solution with ~ Sp ^ Q (i-e., a ray that stays in the equatorial 
plane), but there are also solutions in which ^ and Sp are nonzero. 

Before evaluating the quasi-equatorial light bending, we need to relate the light ray's coordinate distance of closest 
approach, tq, to the invariant impact parameter 6. The distance of closest approach is given by the solution of f = 0. 
From eqn. (jBSP this is a simple quadratic equation in 6, whose positive real solution is 



-= (G- 



2m, 



-1/2 



Alternatively, f = is a cubic equation in rg, whose one real solution is given by 



6 ~ 

Taylor expanding in m,/6 ^ 1 yields 



2 

3V2 



Gi/2 cos 



1 



_33/2 



m. 



G3/2 6 



(Bll) 



(B12) 



6 



-.1/2 





SF-* 




2 


(t) 


2G5/2 


(t) 


G4 



/m,y 105 F** /rn,\ 

Kb') ~ 8gTi75 V~6~J 



3003 f12 /m. 



16G17/2 



(B13) 



(We could further expand F and G as Taylor series in m,/6, but choose not to do that yet.) Note that in the absence 
of spin (a = 0), F = G = 1 and so eqns. (jBlll) - (jB13l) reduce to their respective Schwarzschild values in [17]. 



2. "Horizontal" Light Bending Angle 

We consider the bending of a null geodesic along the (/)-direction (horizontal). From eqns. (jBSp and (IB9[) . we can 
write the equation of motion as 

^^i^, s6ri/2 (r-2m,F) 

dr f [a2-f r(r-2m,)][r3 + 62(2m,F2-Gr)]i/2 ' 
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To understand the sign, consider Figs. [T]and|3l In the case of retrograde motion, 0/ = — i? and 0i = tt + ?95, with 
(f)i> (f)f (cf. Section HIT)) : recall from Section Hill that i? and 1)5 are the respective projections onto the zy-plane of the 
ang les d and ds- For the "incoming" ray segment (from the source the point of closest approach), we have (see, e.g., 
[H p. 189]), 



dr , 



where 0o is the value of at the point of closest approach. For the "outgoing" segment (from the observer the point 
of closest approach), we have 



/•OO 


d0 


/ ro 


dr 



Putting them together yields 



- Cj)f 



7T + ds+'d = (pi-(l)f 



/■OO 




Jro 


dr 




-'J 



dr . 



poo 


d(t> 


' ro 


dr 



dr . 



Identifying i) + as the "horizontal" bending angle dhor (for quasi-equatorial lensing; see eqn. ([H])), we can rewrite 
this equation in the more familiar form (cf. [17|) 



In the case of prograde motion, we have (j)f = d and ipi — — (tt + 1^5) with (j)f > (pi. Similar logic then yields 



poo 




J ro 


dr 







dr 



= 2 



POO 


d(t> 


1 ro 


dr 



dr . 



(B15) 



(B16) 



Identifying d + i)s — dhor again yields eqn. (jBlSp . 

Thus, eqn. (jB15l) represents the general expression for the "horizontal" component of the bending angle. The 
integrand depends on the invariant impact parameter, &, but the integral itself also depends on the coordinate 
distance of closest approach, tq. For pedagogical purposes, and to connect with previous studies of lensing by Kerr 
black holes, it is useful to express the integral purely in terms of ro, and later to convert back to h. 

In the weak-deflection regime, r — 2m, and r — 2 m, F are always positive, so all factors in eqn. (|B14p are positive 
except for s = ±1. Hence the absolute value in eqn. (jB15l) simply removes the factor of s. Changing integration 
variables to x — r^/r, we can write the bending angle as 



dhor = 2 



l-2Vhx 



{l-2hx + a?h^x'^) [G(l-a;2)-2F2/i(l-a;3)]i/2 



dx 



where h = m,/ro, and we have used eqn. (jBllI) to substitute for h in terms of rp. In the weak-deflection regime <C 1, 
so we can expand the integrand as a Taylor series in h and then integrate term by term to obtain 



where 



ahor(/l) 



Co TT 



- 4ci/i + 1^-4 C2 
'130 C4 



64 



IStt 



-o{hy 



122 IStt , 
— C3 - — ds I /i^ 



(B17) 



Co 

Cl 
C2 



Gi/2 



- 1, 

G - FG 



G3/2 

f2(f2 + G-FG) 



G5/2 
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d2 = [l5F4-4G(F-l)(3F2 + 2G)-2a2G2] 



C3 



^ [61 - G (F - 1) (45 F"* + 32 F^ G + 16 G^) - 4 G^ (2 F^ + 2 G - F G)] , 



61 



da = — d2 , 

p2 

C4 = — [65fS-49(F- l)F^G-8F2(-4 + a2 + 4F)G2 + 4(4 + a2(F-2) -4f) G^l , 

65 G-^/^ 

d4 = 5-— [ll55F*- 840(F-l)F^G-140F''(-4 + a2+4F)G2 + 80(4 + a^(F-2)-4F)F2G3 

1155 G ' 

+ 8(16- 12a2 + a* + 8(a2 - 2)F)G^] . 

(Terms beyond order four in the bending angle series can be derived but are not used in our study.) In the absence 
of spin, we have a = and F = G = 1, so the coefficients become 

Co =0, ci = C2 C3 = d2 = da = C4 = d4 = 1. 

In this limit, eqn. (|B17p reduces to the correct Schwarzschild result in ("l7j . 

Let us briefly consider the bending angle to lowest order in m./rp and a/rg. At first order, 6 ~ tq so from eqn. (|Blip 
we have F w 1 — s a/rg and G « 1. This yields Cq « and Ci w 1 — s a/rg. So to lowest order eqn. (|B17[) gives 

ahor ~ 4 — 1 - s — 
ro \ ro 

which recovers the known result for such a regime (see., e.g., 0, p. 281]). 

The expression (jB17p is coordinate-dependent because it involves the coordinate distance of closest approach, tq. 
We must rewrite the formula in terms of the impact parameter b to obtain an invariant result. We use eqn. (jB12p to 
write h — m,/ro as a Taylor series in m,/&. 



-m.\ f2 /m.\2 5f4 f^')^ 

gV^ y~) ^ G2 [TJ + 2GV^ \T) ^ "g^ 
231 F« /m.\5 112 f1" /m.\6 



m 



■ S0.3. (t) - ^ (x)^ °(t) • 

and insert this into eqn. (jB17l) to obtain a series expansion for the bending angle in m,/6: 



(B19) 



where 



Cn = ( ^ - 1 1 TT . 



Ci = 4 



F^ + G-FG 



G' 



2 



C2 = [l5F4-4G(F-l)(3F2 + 2G)-2a2G2] , 

C3 = ^ [i6F6-4G(F-l)(3F^ + 2F2G + G^)-a2G2(2F2 + 2G-FG)] 



C4 = TTTTT [ll55F8-840(F-l)FSG-140F'*(-4 + a2+4F)G2 



-.2 

I I I : I : I i-i 1 \ 1-1 — I 1 ri 1 ... - — 1 1 i-> 1 z_l — i— / 1 —1— ij_ i-i 1 1 

64G13/2 

■80(4 + a2(F-2) -4F) F^ G^ + 8 (16 - 12a2 + a'' + 8(0^ - 2) F) G^] 



Eqn. (jB19p holds for values of F and G where the spin is bounded, a? < 1. In other words, when expanding in m,/b 
we really ought to expand in a/6 as well. Formally, we may accomplish this by writing F and G in terms of a and m,/6 
as in eqn. (IBlip . expanding the coefficients Ci as Taylor series in m,/6, and collecting terms to obtain a new series 
expansion for the bending angle. The result is: 
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where 
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107rsa + 4a , 

2857ra2 

- 192 s a H 4 s a' 
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(Recall that cP' < 1.) When there is no spin (a = 0), the coefficients reduce to Ai — 4, A2 — 157r/4, A3 — 128/3, and 
A4 = 34657r/64 and recover the Schwarzschild values in [l7| . 

Note that in our scaled angular variables (|2ip - ([M)) . the "horizontal" bending angle to third order in e is 

. , , 4 157r + 16 s a - I66I1 2 
aho.(e) = ^2 ^ 

256 + 24a2- 607rsa + 64D26';l- 457r6ii+48sa6li+ 246l?- 246io6l2 . .^,,4 
+ "-^^ e^ + Ois) . 

(B20) 



3. "Vertical" Bending Angle 



This section presents new results on the vertical component of the bending angle in quasi-equatorial lensing. From 
eqns. (|B8I) and (IB10[) . the quasi-equatorial light bending in the C-direction is governed by the equation of motion 



where 



dr r 



±i{r) 



1/2 



(r4-62Gr2 + 2m.62 f2 r)i/2 
The equation of motion has solutions of the form 
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5tp 
01/2 



sin [±/(r) -I- p] , 



where 



I{r) = / i{r')dr' , 



(B22) 



and p is a constant of integration. We are interested in the two asymptotic values (r — ^ 00) 

5lp 



c± 



Gi/2 



sin (±/oo + p) , 



where /oo — limr_>.oo li""")- We can eliminate p and relate the two solutions to one another: 
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sm 



-2h 



sm 



31^2 

5tf 
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The asymptotic values C,± must correspond to the initial and final values, Ci and C/, introduced in Appendix]^ but we 
must determine the correspondence. In order to do that, we first examine Ci and C/ more carefully, using eqn. <\A2\ . 
Recall that in the quasi-equatorial regime we have Lp = Lp^ -\- Sip and (ps — (po + + ^Vs? with (/Sq = tt for prograde 
motion and (po = for retrograde motion. Using these relations, eqns. (IA2p and (jA3p both become 



sin^i = sin??5 sa\5(ps, sinC/ = sini? smSip 
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Since we are working to first order in C and S^p, we can write these as 

Q — Sips sm-ds 7 C/ — Sip sin I? . 

Upon considering the spherical case (see below), we recognize that we want to put C_ = Q and C+ = Cy in eqn. (jB23[l . 
This substitution yields 



Sip 
G^/^ sm 



sm 



-2/oo +sin"^ f G^/^sini? 



sin-iJc 



(5(p . 



(B24) 



Notice that the coefficient of Sp depends only on not on p. (In addition to the explicit i? dependence, there is 
implicit dependence through 'ds and loo, which depends on 6 = c?i sin??.) We can therefore define it to be the function 
W{'d), with a factor of sint^s that will prove to be convenient later. 

Before evaluating loo, let us check the case of a spherical lens to make sure our result is reasonable. For a spherical 
lens, a = and F = G = 1, so we have 



b dr 



ri/2[r3 _62(^_2m.)]i/2 



where the last equality is obtained after comparison with the spherical limits of eqns. (|B14p and (jBlsp . Together with 
our choices C_ — Ci and C+ = Cf, cqn. (jB24p then becomes 



Sps = 



Sp 
sin -fff 



sin [— TT — — I? + sin (sin -i?)] = Sp , 



which is consistent with the symmetry. This verifies our choice of signs above. 

We now evaluate the integral (for the general case, not just the spherical limit), in parallel with the analysis in 
Appendix IB 21 From eqns. (|B21|) and (|B22|) we have 



6G1/2 



(r4 -62Gr2 + 2m.62F2r)i/2 
Using eqn. (jBll|) for b and changing integration variables to a; = j'o /r yields 

.1 Ql/2 



dr . 



/o [G(l-.x2) -2F2/l(l-2;3)]l/2 

where h — m./rg. Taylor expanding in h and integrating term by term gives 



dx , 



2F2 



- + h H ^(ISvr - 16) t 

2 G 8G2^ ' 12 G3 



(244 - 457r) + ——(-1664 + 6937r) + O [hf . 



128 G4 



We now use eqn. (|B18p to write h in terms of m,/6, and then collect terms to obtain 



2f2 /m, 

G372 



ISttF^ /m,^ 2 
8G3 VT 



64 F^ /m. 

3G972 



34657rF^ /m. 



128 G6 



m. 



As in eqn. (|B20[) . when we expand in m,/b we ought to expand in a/6 as well. We use eqn. (IBlip to write F and G in 
terms of a and m,/6, and then collect terms to find 



TT 

2 



/ 34657r 
V 128 



15n 



4sa 



128 s a 



1357ra2 
8 



m. 



6 s a' 



64 ISttso 

y 2 



m. 



(B25) 



This is to be used with eqn. (|B24p to describe the "vertical" bending (see eqn. ([^ in Section ITV Al above) . Note also 
that the expression inside the square root in eqn. (j27p is 



l-Wi-df = 1- 



1 



Gi/2 



-2/oo +sin"i (G^/^gij^^ 



(1 -sin^t?) +4sin2i? (^) - 



16 cos 2i!) 



IStt 



i2t? 



m. 



(SOtt - 64 s a) cos 2i9 + a (157r s - 10 a) sin 2i? - 4 tan 
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Since < < 7r/2 and m,/6 <C 1, eqn. (jB26|) is nonnegative. 

Finally, analogously to the "horizontal" component of the bending angle derived in Section IIIIl we derive the 
"vertical" component of the bending angle, as follows. Consider the angles Vi and vj shown in Fig. [TJ We define vj 
to be strictly nonnegative and within the interval [0, 7r/2), but allow Vi to be negative, so that — 7r/2 < Vi < 7r/2, and 
enforce the following sign convention for Vi. As shown in Fig. [1] Vi is the angle whose vertex is the point B' on the 
lens plane, and is measured from a line parallel to the equatorial plane. If Vi goes away from the equatorial plane, 
then we take it to be positive; otherwise it is negative (e.g., the Vi shown in Fig. [T] is positive). Now denote by Vi 
and Vf the respective projections onto the xz-plane of the angles Ui and Vf, and adopt the same sign conventions for 
them. With these conventions, the "vertical" component of the bending angle can be unambiguously expressed as 



ttvort ^ Vf - Vi . 

By the positivity of z>/ and the fact that the bending is nonnegative, we have 

Vi <Vf . 

(Indeed, with our signs conventions the condition Vi > Vf would be equivalent to repulsion of the light ray.) Writing 
Vi and Vf in terms of the angles i^, (p, tps, we have 

Vf = tan~"'^(tan??sin(p) , 

Vi — tan~"'^(tani?s sin(7r — (^s)) , 



which in the quasi-equatorial regime reduce to 



Vi 



=F5(p5 tani?5 



where we have set ip = ipQ + Sip, 5pg = p^ -\- t: + Sps, with pq = (retrograde motion) or tt (prograde motion), and 
expanded to linear order in the small angles dp and dps- Using the identities W{'&) dp = siwds dps and 7?s = cihor — 
given by eqns. (|B24p and ([TT|). we can thus write dvcrt as 



avert ~ ±dp 



tani9 + 



COs(dhor - ??) 



The expression inside the square brackets is of the form 16D cscd sec^i?e^ + 0{e)^, so it is positive (recall that 



< i} < 7r/2). Since the bending angle is strictly nonnegative, we will adopt 
so that we may write 



Svort ~ dp 



tani? 



COs(dhor - 19) 

We now expand eqn. (|B27|) in our scaled angular variables (PT|) - ([M)) to third order in e to obtain 



for dp>0 and " - " for dp < 0, 

(B27) 



Avert (e) 



dp < tan i3 H — 777- sin 
1 Gi/2 



-2/00 +sin"^ (G^/^gjj^^ 



1 



COs(Q!hor - '&) 



, 4 157r-32sa - 166*1 2 

'^1^^^ ' 

256 + 72 a2 - 907r s a + GiD^Of^ + 96 s a 61- AbirOi + 2461^ - 246*0612 



0{e) 



(B28) 



The result in eqn. (IB28I) is new. 
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Appendix C: Quasi-Equatorial Time Delay 



We now compute the time delays for quasi-equatorial lensed images. Let i?src and i?obs be the radial coordinates 
of the source and observer, respectively. From geometry relative to the flat metric of the distant observer, who is 
assumed to be at rest in the Boyer-Lindquist coordinates, we can work out 



i?obs = dL , i?src — {d^s 



4tan2 B)'^\ 



The radial distances are very nearly the same as angular diameter distances since the source and observer are in the 
asymptotically flat region of the spacetime. In the absence of the lens, the spacetime would be flat and the light ray 
would travel along a linear path of length ds/ cos B from the source to the observer. 

The time delay is the diff'erence between the light travel time for the actual ray, and the travel time for the ray the 
light would have taken had the lens been absent. This can be written as 



cr = T{R,,,) + r(i?obs 



with 



We use t and f from eqns. (IB7p and 
X = ro/r. This yields 



T{R) = f 

J rn 





dt 




i 


r 




dr^ f 






dr 


Jra 


r 



ds 
cosB 



dr . 



T(i?) = ro 



, substitute for b using eqn. (jBlip . and change integration variables to 
(G- 2f2 /i)i/2(i + ^2 /i2 2,2) _ 2saF/i2 



,/jiX^{l~2hx + a? x^) [G(l - a;2) - 2 f2 h{l - x^)]^' 



dx , 



where h — m,/ro. We expand the integrand as a Taylor series in h and integrate term by term. The result is a series 
in h whose coefficients are rational functions of a; = ro/ R. The first three terms in the expansion are 



T{R) = ^R-'-rl + hro 



G(l + 

3f4 + 4f2g + 8G2 



2 In 



(Cl) 



2G2 



sin ^ (jj 



0{hY 



f2 (4 G + (f2 + 4 G)w) Vl -cj2 
2G2(1 + cj)2 

The third-order term is easily obtained, and is needed in the derivation of eqns. (|48p and (jSOp . but is too unweildy to 
write here. Note that if we substitute for F and G using eqn. (|Blip and take the far-field limit, we recover previous 
results (e.g., [Hill). 



Appendix D: Third-Order Results 



The vanishing of the third-order term in the "horizontal" lens equation (|25p yields 

^3 = ^aoa2]n2 , , J 128 s a [-384 + 9(57r - 2gi)gi + 4^0 (3^2 + 8go(2iJ(3 - 2 jJgg) 
7b8fc'o (t'o + 1) [ 

+ + D){-1 + D6l)S))\ + 768a^(27r-30i)- 768sa3 

+ 157r [3(487 + 4861^) + 326io (-3^2 + 86io(2L>(-3 + 2Del) - 3(-l + D){-1 + Del)£))] 
+ 2566*1 [-3(48 + Oj) + 26*0 (36*2 + 49o{D{6 + DOliJ - 661)) - 3(-l -I- D){1 + 2D0l)^))] 
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The third-order magnification term in eqn. (|34p is 

1 



M3 



12288 [el - 1)4(00 + 1)^ 



a Ml + Ms + M3 + M4 



where 



Ml = 32s6lo 



+ 20AS{el - l)^(-3 + 3(-l + D)Del^£, + C (2(-3 + D)D 
+ 3(-l + i:')(l + 4D) + 6'o(-3(4 + + D{-Q + 2D + 3^)) 
+ 26i6(-6 - 9e - 97:>(2 + + 2I?2(8 + 9^) + fo(33 + 6^ - 6L'(-2 + £'(4 + C))) 

+ 30^(2(3 + ^) + D{A + 5e - i?(8 + 70)) 

M2 = 3847r 6*0 [(6'o - 1) (15 + 956*^ + 2316ifS- 179761^5 + 17776ig- 1287C + 22176*^^ + 296'^'*)] , 
M3 = -49152s6'^(5-66';^ + 186';5-60^ + 5e'g) , 

M4 = 457r6lo(6lo -1){-16(-1 + 61^)2 [-25 - 128L>2 + (-125 + 128(4 - 7i:')L>)6l^ 

+ (-1711 + 128L'(-8 + 5D))e^ + (1461 + 128(8 - l5D)D)e^ + 256(-2 + D)De^] 

- 2257r2(-l + 02)2(1 + 7(91 + 39^ + 59^)) 

- 4096(-l + D){9o - 9l)^{-l + 9l{2 - 29^ + 0^ + d{9 + 9^ + 50^ + 06)))^} . 

(We will forgo writing as the expressions are too lengthy.) Expressed in terms of /3, the third-order relations are 
+ ^ _ _ s+ a 2257r2 - 2048(1+ (-1 + £»)!) 



128 |/3| 

TT a2 (-7168 - 1224;32 ^ 1272^* 230;36 + 15/3^) 157r M5 



32 /32 (^2 ^. 4)7/2 4096 (/32 + 4)V2 

+ _ _ s+a Ms 57ra2 (3/?^ + 4) ^ ,3 16 + 14 + 5/3^ 

M3 - /^3 - 384 /32 (/32+ 4)7/2 + p 8s a ^4 (^2 + 4)7/2 

157r (25(-16 + 97r2) + 2048L»(-L» + 2(-l + D)^) 
40961^1 ' 



(Dl) 



where 



Ms = 22577^(140 + 70/32 ^ 14^4 _^ ^6) _ i6(-4 ^ ^2) (3g72 + 250;32 25/3'' - 512(1 + /32) ^) 
+ 4096i:>(4 + /32) (4 + lOe + /3^(4 + (8 + /32)e)) 

- 20481)2(4 + /32) (12 + 24^ + P^{U + jS^ + 2(10 + /32) ^)) , 
Ms = -61447:>(4 + /32) (8(2 + ^+(3^(2 + (3^) (4 + (6 + P^) O) 

+ 20481)2(4 + /32) (-24 + 48^ + /32(4(-8 + ^2) + 3(2 + /32)(8 + /32) f)) 
+ 3 [2257r2(8 + 92/32 ^ 70^4 _^ ^4^6 ^ ^8) 

- 20 48(4 + /32)(/3*^ + 8(-2 + 0+ 4/3^(2 + C) + 2/3^(5 + 0)] ■ 
Finally, the third-order centroid term is 

aCj + a2 C3 + C4 + C5 



e 



cent, 3 



12288 /32 (/32 + 2)4(/32 +4)5/2 ' 
where 

C2 = 16s;32(2 + /32^(4gQ^(2 + ^2^2(4^ /32)(14 + 5/32) + 



+ 22577274+732(1168 + 2OI6/32 + 1642/3"* + 597/3"^ + 97;3^ + 6/3^°) 

- 2048(4 + /32)3/2(6(/38 + 16(4 + + /3'^(12 + + /3*(46 + 5^) + /32(68 + 6^)) 
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- 6Li(-32 + /3*^(-2 + + 48^ + 5l3\~2 + 30 + 2p^{-6 + 29^) 
+ i:i2(-5/36 + 60/34e + 32(-5 + 60 +4/32(17+780))) , 

C3 = 384P{2 + P^){-S2P^{2 + P^f{'i+ ;32)3/2 

+ 7r(-18176 - 15024/3^ + 1512^^ + 5172/3^ + 1926/3^ + 280^^° + IS;^^^)) , 
C4 = -24576s^V4 + ^2(-32 + 84^^ ^ 44^4 _^ 5^6) ^ 

C5 = -/3^(2 + /32)2(108007rV4 + /32(8 + 6^^ + + 101257r^(268 + 280/3^ + 105/3^ + 17/3^ + ^^) 

- 32768(2 + /32)(4 + p'^f^^{-D{-6 + D{8 + /32)) + 3(-l + £>)(-! + 2Li(2 + p^)) 

+ 2407r(4 + /32)(-3(4472 + 3198,52 + 556/3^ + 25/3^) + 768D(12 + 2/3" + /32(10 - 30 - 100 
+ 768(6 + 5/32 + /34)^ _ I28£)2(i88 + 3^6 _ 34^ _,_ 6/34(7 + ^) + 4/3^(40 + 30))) • 
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